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flame  Is  a  homogeneous  gaseous  atmosphere  with  plane-parallel  stratification, 
(b)  the  gas  consists  of  Inert  molecules  plus  sodium  atoms  which  can  be 

excited  to  the  BP^  or  BP^  level,  (c)  there  is  local  thermodynamic 

Iqul librium  governed  by  the  local  temperature,  (d)  the  temperature  gradient 
can  be  represented  by  a  parabola  whose  vertex  is  at  the  center  of  the 
flame,  (e)  the  dispersion  profile  and  number  density  of  sodium  atoms  have 
average  values,  Inside  the  flame,  that  are  Independent  of  depth,  and  (f) 
the  Individual  line  dispersion  profile  Is  replaced  with  a  two-line  function 
to  simultaneously  describe  the  spectral  distribution  of  both  of  the  sodium 
D  lines. 

The  parameters  of  the  radiative  transfer  theory  were  supplied  from 
calculated  thermodynamic  properties  of  the  flare.  Optical  thickness  as 
a  function  of  position  In  the  flame  was  determined  using  computed  sodium 
atom  densities  and  physical  flame  size  obtained  pootographically.  A 
flame  temperature  gradient  was  constructed  numerically  as  a  function  of 
temperature  In  the  flame  using  the  computed  temperature  at  the  flame 
center  and  the  boundary.  The  two-line  dispersion  profile  was  constructed 
as  a  function  of  line  broadening.  The  shape  and  Intensity  of  the 
broadened  flare  spectrum  was  computed  without  introducing  further 
assumptions. 
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FOREWORD 


This  paper  was  prepared  for  presentation  at  the  annual 
NAVAIR  review  of  Its  sponsored  research  programs  in  the  field 
of  Energy  Conversion  (non-propulsive  aspects)  to  be  held  on 
26-27  March  1974,  at  the  Denver  Research  Institute,  University 
of  Denver.  The  work  Is  sponsored  by  the  Research  and  Technology 
Group  of  the  Naval  Air  Systems  Command.  Dr.  Hyman  Rosenwasser 
is  the  Program  Manager  in  NAVAIR. 

The  progress  report  for  fiscal  year  1972  is  contained  in 
RDTN  No.  223,  "Pyrotechnic  Flare  Spectroscopy  III",  Naval 
Ammunition  Depot,  Crane,  Indiana,  1  November  1972. 
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INTRODUCTION 


Illuminating  flares  are  typically  made  from  a  mixture  of 
magnesium,  sodium  nitrate,  and  a  binder.  Light  is  emitted  from 
these  flares  at  a  luminous  efficiency  of  about  50,000  candle- 
seconds/gram.  To  satisfy  the  continuing  need  to  generate  light 
more  efficiently,  the  specific  objective  of  this  research  was 
directed  toward  determining  the  mechanisms  by  which  light  is 
emitted  from  illuminating  flames,  the  new  knowledge  providing 
the  basis  for  future  improvements.  In  this  paper,  I  will  give 
a  suimary  of  the  recently  completed  work  during  which  we 
developed  and  validated  a  radiative  transfer  model  of  a  pyro¬ 
technic  flare.  The  model  allows  us  to  make  predictions  of  the 
flux  emitted  from  the  flame.  Further  details  concerning  this 
research  can  be  found  in  the  final  report.1 

In  this  paper,  I  will  first  present  the  model  and  its 
theoretical  basis,  next  experimental  data  needed  for  the  model 
validation,  and  finally  a  comparison  of  experiment  with  the 
model . 


THEORETICAL 

It  was  hypothesized  that  the  flux  radiated  from  a  flame  of 
burning  magnesium,  sodium  nitrate,  and  binder  consists  mainly 
of  photons  from  the  resonance  transition  of  the  sodium  D  lines. 
If  this  were  the  case,  we  should  be  able  to  predict  radiative 
output  by  treating  only  the  sodium  D  radiation.  It  is  on  this 
basis  that  the  radiative  model  was  developed. 

We  know  that  the  total  radiant  intensity  Ivy(T)  at 

frequency  v  in  a  direction  described  by  u  =  cosO  and  issuing 
from  a  volume  element  at  optical  depth  x  is  given  by  the 
radiative-transfer  equation3 


dIvu(l)/dT  ■  ♦va[IvtJ(T)  -  S  v(  T )  ]  .  (1) 

The  optical  depth  x  is  related  to  the  physical  depth  z  by 
x  *  fk  dz  ,  where  k  is  the  linear  absorption  coefficient  of 

V  V  v 

the  flame.  The  normalized  spectral  two-line  profile  of  the 
absorption  coefficient  4>va  is  a  function  which  takes  account 


1 


of  flame  line  broadening  mechanisms.  The  line-source  function 
S  (t)  accounts  for  increments  or  decrements  in  the  radiant 

V 

Intensity  from  a  volume  element  at  optical  depth  T  due  to 
emitters  and  absorbers  within  that  volume  element.  It  is 
defined  at  a  given  frequency  by  S  =  E  /k  ,  where  c  is 

V  V  V  V 

the  monochromatic  volume  emission  coefficient. 

Formal  integration  of  the  transfer  equation  yields  the 
expression 
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‘v.  eXpC'(VT,)*wi/,,:i 


r  t =t 


)4>va/ y 3  exp[-(T-x  )<i>va/y ]dx  * 
(2) 


where  t  and  t  are  the  optical  depth  integration  limits 
1  2 

from  front  to  the  rear  of  the  atmosphere  respectively,  and 
Ivi  and  1^  are  the  spectral  intensity  at  optical  depths  T 

and  T  respectively.  In  order  to  solve  the  transfer  equation, 

2 

for  the  observed  radiant  intensity,  the  flame  is  represented  by 
the  following  model. 

(1)  The  flame  is  a  homogeneous  gaseous  atmosphere  with  plane- 
parallel  stratification. 

(2)  The  gas  consists  of  inert  molecules  plus  sodium  atoms 
which  can  be  excited  to  the  aP,  or  ZP,  level. 

£■  -8- 

(3)  There  is  local  thermodynamic  equilibrium  LTE  governed 
by  the  local  temperature. 

(4)  Energy  exchange  by  radiation  leads  to  radiative  equilibrium. 

(5)  The  refractive  index  of  the  medium  is  unity. 

(6)  The  radiation  is  unpolarized  when  emitted  and  remains 
unpolarized  in  its  interactions  with  flame  species. 

(7)  The  temperature  gradient  can  be  represented  by  a  parabola 
whose  vertex  is  at  the  center  of  the  flame. 

(8)  The  absorption  profile  <j>  and  number  density  of  sodium 
atoms  Nq  have  average  values,  inside  the  flame,  that 
are  independent  of  x  . 
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The  form  of  Eq.  (2)  has  been  simplified  for  the  present  case, 
namely  (a)  the  observed  flux  Is  that  emerging  normal  to  the 
surface  (w=l),  (b)  no  flux  is  Incident  on  the  rear  surface  of 
the  atmosphere  ( I vS3=0) ,  and  (c)  Sv(t)  =  Bv(T' )  for  the  LTE  case 

where  Bv(T')  Is  the  Planck  function.  Under  these  conditions, 

integrating  from  the  front  surface,  where  z  and  t  are  0  , 

to  the  rear  surface  where  t  «  T  ,  the  total  optical  thickness, 

the  monochromatic  emergent  intensity  is 

V-O  Bu(T')  exp(-T*va)dr  .  (3) 

J  T  =  0 


The  theoretical  relative  spectral  radiant  power  , 

proportional  to  spectral  emergent  intensity  I  °  for  a  particular 

model,  was  found  by  numerical  Integration  of  Eq.  (3)  using 
Simpson's  rule  of  2m  intervals.  Each  computed  spectrum  $x  , 

normalized  to  a  power  maximum  of  unity,  is  plotted  for  comparison 
with  the  corresponding  experimental  flare  spectrum  in 

Figs.  1  and  2. 

Parameters  of  the  theory  that  must  be  supplied  from  properties 
of  the  flame  are  (1)  optical  thickness  t(z)  as  a  function  of 
position  in  the  flame,  (2)  a  flame  temperature  gradient  T'(z) 
as  a  function  of  position  in  the  flame,  and  (3)  the  scattering 
profile  a  parameter  in  $  . 


EXPERIMENTAL 

The  test  flares  were  composed  of  50  g  of  a  magnesium-sodium 
nitrate-binder  mixture  compressed  into  3.3  cm  i.d.  by  5.5  cm 
long  paper  tubes,  having  formulas  shown  in  Table  1.  Formula 
groups  1,  2,  and  3  are  nearly  stoichiometric  mixtures,  the 
sodium  nitrate  in  groups  2  and  3  being  .1  and  .01  of  group  1 
respectively.  Stoichiometry  was  maintained  in  groups  2  and  3 
by  addition  of  potassium  nitrate  chosen  because  it  reacts  with 
magnesium  at  about  the  same  rate  as  sodium  nitrate  and  because 
of  its  low-emissivity  in  the  neighborhood  of  the  sodium  D  lines, 
the  region  of  Interest  for  these  studies. 
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Each  of  three  dlfferen  Illuminating  composition  formulas 
was  tested  in  air  at  8  levels  of  pressure;  namely  760,  630, 

300,  225,  150,  75,  30  and  6  torr.  For  each  pressure-formula 
combination,  the  burning  time,  flame  size,  and  relative  spectral 
radiant  power  distribution  in  the  visible  region  were  recorded. 

Relative  radiant  power  spectra  4>  '  of  typical  flares  for 
each  pressure-formula  combination  areAplotted  In  Figs.  1  and  2. 
The  solid  curves  in  Figs.  1  and  2  are  the  experimental  data. 
These  spectra  are  normalized  to  a  peak  value  of  unity  for 
convenience  In  the  first  step  of  the  theoretical  comparison. 
Spectra  were  not  obtained  for  formula  groups  2  and  3  at  6  torr 
because  the  flares  did  not  sustain  combustion  at  this  pressure. 
Group  1  flares  at  6  torr  barely  burned.  Combustion  difficulty 
was  visually  observable  for  all  flares  tested  at  75  torr  or 
less. 

Values  of  the  physical  flame  depth  z'  range  from  6  cm  for 
formula  group  1,  760  torr  to  2.5  cm  for  formula  group  3, 
torr,  a  rather  narrow  range  considering  the  large  range  of 
experimental  conditions. 


DISCUSSION 

Visual  comparison  of  the  theoretical  4>x  and  experimental 

4>x'  relative  radiant  power  spectra,  plotted  in  Figs.  1  and  2, 

shows  that  the  distribution  computed  from  theory  agrees  quite 
well  with  the  experimental  distribution  for  each  pressure- 
formula  combination. 

In  summary,  it  has  been  shown  that  the  spectral  radiant 
power  distribution  of  a  pyrotechnic  Illuminating  flare  flame 
can  be  predicted  by  a  two-line  radiative  transfer  model  which 
has  been  described.  This  can  be  done  without  Introducing 
assumptions  which  require  ad  hoo  modifications  of  the  model  to 
describe  different  flares.  Known  system  variables  such  as 
flare  formula,  flare  size,  and  ambient  pressure  are  the  necessary 
and  sufficient  input  needed  for  the  theoretical  prediction. 
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TABLE  1.  Flare  Formulations 


Ingredients 

Formula  Groups 

1  2  3 

Magnesium 

44. 0a 

40.4 

40.04 

Sodium  nitrate 

51.5 

5.15 

0.515 

Potassium  nitrate 

— 

49.95 

54.945 

Epoxy  binder  mix 

4.5 

4.5 

4.5 

aPercent  by  weight 


GROUP  1  GROUP  2  GROUP  3 


Fig.  1.  Illuminating  Flare  Flame  Spectra  for  formula  groups  1,  2,  and 
3  at  4  levels  of  ambient  pressure.  Theoretical  relative  radiant  power 
values  $  are  indicated  by  boxes  (").  Experimentally  determined 
relative  radiant  power  values  are  shown  by  the  solid  line. 
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Fig.  2.  Illuminating  Flare  Flame  Spectra  for  formula  groups  1,  2,  and 
3  at  4  levels  of  ambient  pressure.  Theoretical  relative  radiant  power 
values  4>  are  indicated  by  boxes  (»).  Experimentally  determined 
relative  radiant  power  values  '  are  shown  by  the  solid  line. 


